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Abstract 

The cellular automaton model is used to simulate diffusion and aggregation with disso¬ 
ciation of point particles in 2D. A continuous phase transition is found that separates 
creation of compact aggregates and fractal ones. The transition is the function of pair- 
interaction energy (K;,), type of neighborhood and temperature T. Manifestations of 
the transition in real physical systems are discussed. 


1 Introduction 

In the modern theory of fractal growth the signihcant efforts were dedicated to defining models 
([2], [3]) and developing the theoretical concepts necessary to understand these model and to 
compute their properties analytically [1]. The crucial elements of the simulation models and 
theoretical schemes are: 

- irreversibility; 

- screening and freezing; 

- spontaneous evolution to a fractal structure. 

Moreover, the fractal properties are presumed to be well defined only asimptotically with 
respect to the large time and scale limits. It means that the problem of fractal growth consists 
in considering regions of the system that are very far from growing interface and that will 
not continue to be modihed by further growth. But there are many experimental evidents of 
surface aggregates of different morphologies: fractal and dendritic [4], compact and fractal 
[5], convex and concave ones [6], which were determined by kinetics of 2D growth processes. 
Recent simulations brought to light dependence of aggregate morphology from coverage [7] 
in irreversible diffusion-limited aggregation (DLA). These observations of transition from 
fractal to compact aggregates can be explained by deposition of atoms in internal regions of 
the aggregates. But the passage from irreversible to reversible aggregation was also achieved 
by tuning a value of pair bond energy [8]. It points out an influence of dissociation processes 
on morphology of aggregates. 

Investigation of morphology as a function of deposition conditions during DLA with 
dissociation (DEAD) is of great interest in the spirit of the great importance of building 
low-dimensional nanostructures with desirable geometry. That is why in the following the 
cellular automaton (GA) model of metastable fractal aggregates created by DEAD is consid¬ 
ered (Sec.2). In Sec.3 the dependence of fractal dimensions as a function of kinetics of DEAD 
is presented and the results are discussed in Sec.4. 


2 Cellular automaton model of diffusion-limited aggregation with 
dissociation 

The rules of our CA model are similar to the classic DLA model [2] with the following addition: 
particles are randomly generated beyond a growing aggregate with a constant rate F and 
diffuse with a diffusion rate D ~ exp(—where Em is a migration energy. And they 
can dissociate from an aggregate with the following rate R T)eexp(-where 
El, is a pair bond energy, De ~ exp{—Ee/kBT) is a diffusion rate of particles in aggregate 
(for example, in the simplest case D = De [8]). Local interaction with the nearest neighbors 
is taken into account by summation. Our simulations involved a square lattice with sizes 
ranging from (100 x 100) to (600 x 600) with periodic boundary conditions. Simulations 
were conducted for different values of the ratio n = exp{{Ei, + Em — Ee)/kBT) and the ratio 
D/F ranging from 10^ to 10®. The CA evolves in discrete time steps and the state of each 
site at the next time step is determined from the state of the site itself and those of the 
nearest-neighbor sites (NNs). The von Neumann neighborhood (with four NNs) and Moore 
one (with eight NNs and diagonal coupling) were used. 

3 Metastable fractal aggregates with quasi-stable fractal 
dimensions 
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Figure 1: Typical morphologies created by diffusion-limited aggregation with dissociation 
(DLAD) on initial stages for the von Neumann neighborhood at D/F = 10® for (100 x 100) 
sections: a) k = 20; b) k = 100; c) k = 200; d) k = 1000. 



Fig.l illustrates typical morphologies for the von Neumann neighborhood on initial stages 
at D/F = 10® for (100 x 100) sections: a) k, = 20; h) k = 100; c) k = 200; d) k = 1000. 
Fractal dimensions were determined from hts to the linear region in log-log plots, which were 
obtained by density-density correlation function (de) and box counting measurement {db). 
For this purpose a central part of an aggregate was considered as an inscribed square of 
maximum dimension. On initial stages of aggregation (the number of particle in aggregate 





Table 1: Fractal dimensions {dc) for the von Neumann neighborhood obtained by density- 
density correlation measurements. Mean values of the interpolation lines are denoted by word 
’’mean”, standard deviations are shown in row ”sd” and correlation coefficients are given in 
row ”corr”. 


K 

20 

40 

100 

150 

200 

400 

1000 

2000 

mean 

1.956 

1.869 

1.857 

1.794 

1.755 

1.669 

1.653 

1.652 


0.012 

0.014 

0.005 

0.008 

0.003 

0.007 

0.004 

0.006 


0.9998 

0.9997 

0.9999 

0.9999 

0.9999 

0.9996 

0.9999 

0.9996 


Table 2: Fractal dimensions {dc) for the Moore neighborhood obtained by density-density cor¬ 
relation measurements. Mean values of the interpolation lines are denoted by word ’’mean”, 
standard deviations are shown in row ”sd” and correlation coefficients are given in row ”corr”. 


K 

4 

10 

20 

40 

100 

200 

400 

mean 

1.945 

1.862 

1.776 

1.71 

1.708 

1.652 

1.662 

sd 

0.007 

0.007 

0.003 

0.007 

0.006 

0.018 

0.027 

corr 

0.9997 

0.9998 

0.9999 

0.9997 

0.9997 

0.9992 

0.9980 


n < 10^) length scale is very small for reliable calculation of fractal dimensions. In the 
range 10^ < n < 10^ the dimensions grow and achieve asymptotic values for n ^ 10^. The 
values of dc shown in Table 1 for the von Neumann neighborhood and Table 2 for the Moore 
neighborhood with the main statistical characteristics of interpolation lines. Box counting 
measurements demonstrate several linear regions in log-log plots: from initial region with 
df, = 2 to a staircase part with alternating regions of d;, ~ 1 and db ^ 2. The regions have 
small length scales for reliable determination of df,. Although a core of aggregates becomes 
compact with time, they continue to generate ”DLA-hke” branches. When deposition was 
ceased, aggregates underwent signihcant restructions. The previous simulations shown that 
aggregate splitted in many small parts and coalesced with creation of one large compact 
aggregate with dc = 2. 

4 Simulated phase diagrams as a new tool for investigation of 
surfaces 

In the well-known work of Witten and Sander [2] DLA with sticking probability, which de¬ 
pends on the number of bonds to hlled neighbor sites B, was considered without dissociation. 
Particle can stick to an aggregate only during a fraction s of its encounters: s = t^~^, where 
t < 1. In addition to their conclusion that ’’long-range properties seem to be unaltered” for 
0.1 < t < 0.5 we should like to emphasize the crucial changes in DLA caused by dissociation. 
It should be noted that the whole aggregate can split in unlinked parts, but scale invariance 
dehned by density-density correlations remains. 

The aforementioned results allow to propose the new method for investigation of adatom- 
adatom and adatom-surface interactions by means of investigation of scale invariant aggre¬ 
gates. For this purpose careful simulations of DLAD for different deposition rate, k, type of 
interaction are necessary. It will give model phase diagrams for any kind of real aggregating 




































































adatoms and substrate. Testing the real aggregation processes for different deposition con¬ 
ditions and comparing theirs with the simulated model phase diagrams will allow to bring 
to light new features of local interaction at surfaces. We used the naive Arrhenius-type ex- 
pessions for dissociation rates. It is obvious the more realistic scheme will need the more 
sophisticated algorythms and codes with more reliable results. 

In addition to this we recently observed the very complicated band patterns in metals dur¬ 
ing plastic deformation of metals [9]. They can be explained by development of disturbances 
of concentration of point-type defects and their agglomerates [10]. The crucial point of this 
assumption is based on possibility of point defects to agggregate in scale invariant structures 
with different dimensions and consequently different abilities to capture point-type defects. 
Fractal dimensions of aggregates are dependent on pair interaction energy and they dehne 
the selection rate among heterogeneities with further band pattern formation. The extended 
version of aforementioned investigation in 3D will also allow to make clear mechanisms of 
pattern formation in plastically deformed metals. 

Acknowledgments. The authors gratefully acknowledge Prof.A.A.Snarski for helpful com¬ 
ments on this work. 
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